Transport calculations in complex materials: A comparison of the Kubo formula, the 
Kubo-Greenwood formula and the microscopic response method 
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Recently we have introduced the microscopic response method (MRM) to compute the conduc- 
tivity and Hall mobility for complex system with topological and thermal disorder, which is more 
convenient than the Kubo formula. We prove that for a canonical ensemble the MRM leads to the 
same expression as the Kubo formula. When the gradient of carrier density is small, the MRM 
reduces to the widely used Kubo-Greenwood formula. 
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Linear response theory [1 is a rigorous and compli- 
cated procedure used to compute transport coefficients. 
It constructs the observable macroscopic response by av- 
eraging the operator of microscopic response over the 
density matrix of system. Then, to obtain transport co- 
efficients, one has to calculate an imaginary time inte- 
gral resulting from the commutator between the density 
matrix and microscopic response |3]. If a system has sev- 
eral types of elementary excitations and the interactions 
among elementary excitations are strong, it is difficult 
to fulfil the imaginary time integral with controllable ap- 
proximations. For example, the imaginary time integral 
prevents researchers from computing all the important 
contributions in conductivity and Hall mobility of small 
polaronsjl-Q. Amorphous semiconductors and semicon- 
ducting polymers require approximations beyond small 
polarons: the low-lying excited states often contain both 
localized states and extended states [TJj, and the electron- 
phonon interactions in localized states are much stronger 
than those in extended states [lllll2j . Even for the lowest 
order self-consistent approximationjjljj, it is difficult to 
apply the Kubo formula, and consistently include all im- 
portant contributions for conductivity or Hall mobility. 

To compute the transport coefficient for a mechanical 
perturbation, the microscopic response method (MRM) 
is more convenient fl3l [Tij than Kubo formula. A me- 
chanical perturbation such as the coupling with an ex- 
ternal field can be expressed via additional terms in the 
Hamiltonianp], 0|, and the wave function \P' of system 
in an external held at a later moment is determined 
by its initial value and the time-dependent Schrodinger 
equation. The microscopic response can be expressed in 
terms of the changes in wave function induced by the 
external field [H, [lj]. The ensemble average and coarse- 
grained average needed to compute the macroscopic re- 
sponse (transport coefficient) can then be carried out at 
the final stage. Thus for a mechanical perturbation, we 
are able to avoid the imaginary time integral in the Kubo 
formula [H, H|- 

The Kubo-Greenwood formula (KGF)[lj| has been im- 
plemented in many ab initio codes to calculate the dielec- 
tric function and AC conductivity. However the KGF is 
based on a simplified expression for the current density, 



which is borrowed from the kinetic theory of gas. 

In this paper, we prove that for mechanical perturba- 
tions, the MRM is equivalent to the Kubo formula. To 
fulfil this aim, we first write out the observable macro- 
scopic current density to first and second order in exter- 
nal field in the Kubo formulation. Then the same proce- 
dures are carried out with the MRM[14]. We see that 
the macroscopic response calculated in the two methods 
are the same. We discuss the connection between den- 
sity matrix and transition amplitudes at different order of 
perturbation. We will see why the MRM allows easy clas- 
sification of transport processes and expression for trans- 
port coefficient compared to the Kubo formulation[l|. 
We show that the current density used by Greenwood 
is justified only when the gradient of carrier density is 
small. 

In this work, we use the Schrodinger picture. Con- 
sider a system with N electrons and Af nuclei in an 
electromagnetic field with potentials (A,</>), at time 
t, the many-electron state of system is described by 
^'(ri, T2, ■ ■ ■ ,TN;t). To save space, we will not write 
out the nuclear coordinates explicitly. >F' satisfies the 
Schrodinger equation 



ihdV'/dt = H'V, H' = H + V(t), 



(1) 



where V{t) is the interaction between the system and ex- 
ternal field. The time dependence in V(t') comes from 
the external field. H is the Hamiltonian of the sys- 
tem without external field. We use \m) or ^> m and E m 
to denote the m th stationary state and the correspond- 
ing eigenvalue of the N electrons + M nuclei system: 
H\m) = E m \m). If the system is in a thermal bath at 
temperature T, then before introducing V(t), the equi- 
librium density operator is 



\m)P m {m\, P„ 



(2) 



where Z = ^2 n e~@ En is the partition function. In an 
external electromagnetic field, the velocity operator v.; 
for the i th particle is[l6| 



m 



-*[!>, -eAfc;*)], 



(3) 
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where and = —ihV ri are the position and momen- 
tum operators of the i th particle, and e is the charge of 
electron. Because velocity and position cannot be simul- 
taneously measured, one has to symmetrize the velocity 
and position operators in the current density operator. 
One may conjecture that the current density operator at 
point r is [17j[ : 



N 



(4) 



In the MRM[H 0, [3, we avoided J(r). Now we 
show that Eq.flU) leads to a proper microscopic current 
density [l3l. Ilij]. Because a mechanical perturbation can 
be expressed with additional terms in Hamiltonian, the 
states at time t can be described by a wave function 
which is determined by the initial conditions. The mi- 
croscopic current density at time t and point r in state 
*'(ri,r 2 , ■ • • ,r N ;t) is 



j m (r;i) = / dr*'* j(r)tt', 



(5) 



p + p {1) {t) + p (2) (t) + ■■■ . To first order in V(t), the 
deviation p^'(t) from p isQ 

( m | /t ,( 1 )(i)|n) = i C dt 'St-t'){E n -E m) /K (g) 
in ./_ ^ 



(m\V(t')\n)(P n -P m ). 

The conductivity can be read off from the macroscopic 
current density: 

j^(r,t) = Y / Hp {1) (t)\n}(n\%)\m}. (9) 

mn 

To 2nd order in the deviation p^(t) is 

(m\p^(t)\n) = -^e^- E ^Y,{ ( 10 ) 

fe 



t r-t' 

I / 14.11 



OO o —00 



(P„ / dt' / dt" + P m / dt" / dt 



t r t" 



00 j —00 



where dr = dridr', and dr' — dr% ■ ■ ■ dr/y. Integrating 
by parts, one has 



jm(r;t) 



e 2 N 



■m 



L (r) J dr' 



(6) 



Jt"(E k -E n )/h it'(E m -E k )/h 



m\V(t')\k){k\V(t")\n) 



-( / dt' f dt" + [ dt" [ dt') 



OO J — CO 



OO J — CO 



+ / dr'(*'V r *'* - #'*Vr#')> 

2m 7 

where the arguments of W in Eq.© are (r, r2, • • ■ , rjv; t). 
Eq. © has been independently derived from the principle 
of virtual work|18|. the continuity equationfTiJ and the 
polarization der^tyd. The current, operate given m 
Eq.(j4|) is correct, and it will bridge the Kubo formulation 
and the MRM. 

To write out the macroscopic response in the Kubo 
formula, we notice that the time evolution for a system 
involving mixed states is included in the density matrix. 
The basis set should be a group of wave functions without 
any time dependence (l9l| . In the m) representation, the 
matrix elements of j(r) are 

(n|j(r)|m) = — J dr'[# m V r ^ - *;V r * m ] 
Ne 2 f 

— A(r) / dr'^m, (7) 
m J 

where the arguments of ^ m and are (r, r 2 , • • • ,rjv). 
With p in Eq.© as the initial condition, one can use 
perturbation theory to solve the Liouville equation to any 
order in V(t). The density matrix at time t is p'(t) = 



e «"(«*-B»)/V*'( f, »-- E «'Vfc( m |v(t')|jfe)fl l .(jfc|V r (t")|n>}, 

The 2nd order macroscopic response j( 2 )(r, t) is obtained 
from Eq.Q by replacing (m\p^ (t)\n) with (m\p( 2 '(t)\ri). 
We are going to compare Eqs. (|8l9llU|) and ^(r,t) with 
the corresponding quantities in the MRM. 

In the MRM, the macroscopic response is given by [l3|, 
d: 

j(r,t)=^P n (K(*)lJ(r)IK(*))- (11) 

n 

If the initial state is <I> n , then the state ^' n (t) of a system 
at time t in an external field can be determined [19( by 
applying perturbation theory to Eq.([T]): 

*;(*) = a<°>(n,t)* n (12) 
+ V" a (1) (mn, t)^ m + V" a (2) (mn, t)^ m , 

m rn 

where 

a^(n,t)=e- iE ^ h , (13) 

and 

a {1 \mn,t) = ^ l -e- lE ^' h [ dt' (14) 
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e *( B «.- B »)*7»( TO |y(t')|„), 

and 

k k J -°° 

{m\V{t')\k) f dt"e l(Ek ~ E - ]t " /h {k\V{t")\n). (15) 

J — oo 

To first order in V(t), the macroscopic current density 

is 

jW(r,t) =^P„[(*l°)(i)|J(r)|*l 1 )(t)) (16) 



+(^)(t)|j(r)|*( )(i)}]. 

Substituting Eqs. (|12ll3ll4j) into Eq.flTH]), and using the 
fact that V(t') and j(r) are Hermitian operators, one 
finds the same result as Eqs. (18191) . To second order in 
V(t), the macroscopic current density is 

j^M) =^P„[(^°)(t)|j(r)|*( 1 2 )(t)) (17) 



+ (¥^(m(r)\^(t)) + (tf«(i)|j(r)|*«(i))]. 

Substituting Eqs. <|12ll3ll4ll5[> into Eq. ([T7|l . the first 
term of Eq. ([T7l) is the same as the 1st term of j^ 2 ) from 
Eq. flTUl) . the 2nd term of Eq. dTTj) is the same as the 
2nd term of j' 2 ' resulted from Eq. (|T0[) . One can see 
that the 3rd term Eq. (fT7j) equals the sum of the 3rd 
term and the 4th term of j' 2 ' from Eq. (jT0|) . if one no- 
tices: (1) three integrands are the same; (2) the 3rd 
term in Eq. (|17p is a two-dimensional integral in domain 
[— oo, t; — oo, tj; (3) the 3rd term of j^ 2 ) is a successive 

integration fl^dt' J^dt"; (4) the 4th term of j (2) is 

a successive integration J^^dt" dt' . The procedure 
is easy to proceed to any order in held. The equation 
(fTTj) docs not use any specific property of electromag- 
netic field, the procedure works for any mechanical per- 
turbation. Introducing current density operator (j4|) is 
the key for the proof. In the original MRM, one does not 
need current density operator, the macroscopic response 
is obtained by averaging over the microscopic response 
© over canonical distribution. Eqs. (141111) established a 
connection between two methods. 

It is worthwhile to find the connection between the 
probability amplitudes in Eqs. (|131ll5[) and the density 
matrices in Eqs. (181101) . The element of the density matrix 
is the average the product of two probability amplitudes 
over the M members in an ensemble 19] : 

P,nn = j^^2a*Jn,t)a a (m,t), (18) 

a 



where a is the index of a member in the canonical en- 
semble. To the first order in V(t), 

P$n = ^E°»° ) ( n '*) » ) ("*'*) ( 19 ) 
a 

+^E< (1) (™>^ 0) (™> t )> 

a 

where a^ ' (n, f) is the zero order transition amplitude 
from initial state |n) to final state \n), aa\m,t) is the 
first order transition amplitude from initial state \n) to 
final state m), cS®\m\ t) is the zero order transition am- 
plitude from initial state \m) to final state \m), a^'{n,t) 
is the first order transition amplitude from initial state 
|m) to final state \n). With these explanations, substitut- 
ing Eqs. (|13ll4j) into Eq.flT9"]). one reaches Eq.©, which 
was obtained from Liouville equation. To 2nd order in 
V(t), 

P^ = ^EK (0) (M)ai 2 W) (20) 

a 

+a*W (n, t)a<?> ( m , t) + a*« (n, t)a$ (m, t)}. 

In the 1st term of (|20|) . the initial state is \n), a"a\m, t) is 
the 2nd order transition amplitude through intermediate 
states \k). By means of Eqs. (|13ll5[) . the first term of 
Eq. ([20|) is the same as the 1st term of Eq. (fT0|) . In the 

second term of (|20|) . the initial state is |m), a« is 
the 2nd order transition amplitude through intermediate 
states |fe) to final state \n). The second term of Eq. (f2"0")) 
is the same as the second term of Eq. flTUl) . In the 3rd 
term of (|20l) . two final states \n) and \m) come from a 
common initial state \k), all states {\k}} satisfy k ^ n 
and k ^ m can be taken as the initial state. In terms 
of the same trick in comparing the 3rd term in Eq. (|17[) 
and the sum of the 3rd and 4th terms of y- 2 ' (r, t) derived 
from Eq. flTU]) . we can see that the 3rd term of (|2"0"]) is the 
same as the sum of the 3rd term and the 4th term in 
Eq.flTUJ). 

We explain why the MRM is simpler than the Kubo 
formula for mechanical perturbations. To a given order 
in residual interactions, various transport processes con- 
tribute to a specific transport coefficient. In the MRM, 
each transport process is composed of several elemen- 
tary transitions caused by external field and by residual 
interactions [l3l 20]. Because the microscopic response 
is expressed by the wave function of system in external 
field rather than density matrix, each elementary tran- 
sition appears as a transition amplitude (T3|. According 
to Eq.flTU]), t ne state at t is a linear superposition of the 
various order changes induced by the external field. By 
means of Eqs. (|6ll6ll7|) . gradient operator connects two 
components of the final state [13]. In addition, the transi- 
tion amplitude of a higher order transition is constructed 



4 



by first making a product of the sequence of first order 
amplitudes of elementary transitions and then summing 
over all intermediate states. We can depict each trans- 
port process with a diagram, which has one line con- 
necting two components of final state, and several other 
lines presenting elementary transitions. To a given or- 
der of residual interactions, the topology of diagrams 
can help us classify and construct all possible trans- 
port processes [l^. In the Kubo formulation, all time- 
dependence is included in density matrix, cf. Eqs. (|9l7p . 
To a given order in external field, the change in density 
matrix involves different members of the ensemble, cf. 
Eqs. (|19l20p . Besides, the density matrix is bilinear in 
transition amplitude. Therefore for a transport process 
with more than one elementary transitions, one cannot 
express it as a product of propagators. 

Greenwood derived his conductivity expression 
fromfilij 

j(t) = eTr{p\t)v}, (21) 

where p'(t) is the density matrix of system in external 
field. We show that Eq. (|2~Tj) is justified only when the 
fluctuation in the spatial distribution of carriers is small. 
Using Eq.Q and the commutation relation between 
and P i; Eq.Q becomes 

N N 

j(r) = e]T S(r ri )v, - ^ ^[V rs <5(r - r,)]. (22) 

i=l i=l 

Averaging Eq. {22} over state <F' (r i , r 2 , • • • , r N ; t ) , Eq. © 
gives another expression for microscopic current density: 

j(r;i) = ^V r n'(r) (23) 



+eN / dr'^m-^-iWr - eA(r)]f 



where the arguments of in Eq. (l2"3")l are 

(r,r 2 ,--- ,r N ;t), n'(r) = iV/dr'*'*V r *' is the 
number density of electrons at point r in external field 
(A, (f). The 1st term of Eq. ([23|) can be neglected only 
when the gradient of carrier density is small. Using the 
corresponding relation between the MRM and the Kubo 
formulation, the 2nd term of Eq. (|23l) reduces to Eq. (f2"T|) . 
It is obvious that Eq.([2~Tj) and the consequent KGF 
are not suitable to the localized carriers in amorphous 
semiconductors and the d and / electrons in strong 
correlated systems. For a given error Ax of position, the 
error of velocity is Av ~ m~ l h/ Ax. The error of current 
density is n'eAv ~ TO _1 Sn'/Aa; ~ m _1 W r n'. If the 
charge density is uniform like nearly electron gas, the 
uncertainty Aa; of position is infinity, the momentum of 
electron is completely determined Ap x — 0. One can use 
kinetic expression for current density j = 2e^ k Vk/k, 
where /k is the distribution function. 

In summary, we proved that for a mechanical pertur- 
bation the microscopic response method is equivalent to 
and simpler than Kubo formula. To compute transport 
coefficients for mechanical perturbations, the microscopic 
response method is advantageous because of the ease of 
obtaining expression to a given order of residual interac- 
tions consistently. When the gradient of carrier density 
is small, the strict current density Eq. (l23l) reduces to the 
kinetic expression (|2~Tj) . 
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